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ABSTRACT 

The equilibrium structure of homogeneous turbulent shear flow is investigated from a 
theoretical standpoint. Existing turbulence models, in apparent agreement with physical 
and numerical experiments, predict an unbounded exponential time growth of the tur- 
bulent kinetic energy and dissipation rate; only the anisotropy tensor and turbulent time 
scale reach a structural equilibrium. It is shown that if vortex stretching is accounted for 
in the dissipation rate transport equation, then there can exist equilibrium solutions, with 
bounded energy states, where the turbulence production is balanced by its dissipation. 
Illustrative calculations are presented for a k — e model modified to account for vortex 
stretching. The calculations indicate an initial exponential time growth of the turbulent 
kinetic energy and dissipation rate for elapsed times that are as large as those considered 
in any of the previously conducted physical or numerical experiments on homogeneous 
shear flow. However, vortex stretching eventually takes over and forces a production- 
equals-dissipation equilibrium with bounded energy states. The validity of this result 
is further supported by an independent theoretical argument. It is concluded that the 
generally accepted structural equilibrium for homogeneous shear flow with unbounded 
component energies is in need of re-examination. 
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No. NAS1-18605 while the author was in residence at the Institute for Computer Applications in Science 
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1. INTRODUCTION 


Homogeneous turbulent shear flow has been the subject of a variety of experimental, 
computational and theoretical studies during the past four decades. The popularity of 
this flow lies in the fact that it accounts for an important physical effect — the alteration 
of the turbulence structure by shear — in a simplified setting unencumbered by such 
complications as rigid boundaries and mean turbulent diffusion. Von Karman 1 first 
proposed the problem of homogeneous shear flow which gave rise to some mathematical 
studies during the 1950’s (c.f. Townsend 2 and Hinze 3 for a review). It is a difficult 
flow to simulate experimentally and the first tangible suggestion for its generation was 
put forth by Corrsin 4 . The first successful experimental realization of homogeneous 
shear flow in the laboratory was achieved by Rose 5 and was followed by a series of 
landmark experiments by Corrsin and co-workers 6-8 . More recently, Tavoularis and 
Karnik 9 and Rohr et al. 10 performed more exhaustive measurements of homogeneous 
shear flow shedding new light on its basic structure. 

With the dramatic increase in computer capacity achieved by the late 1970 ’s, direct 
numerical simulations of homogeneous turbulent flows became possible. Rogallo 11 was 
the first to conduct a direct simulation of homogeneous shear flow; the results that he 
obtained were well within the range of the experiments of Tavoularis and Corrsin 8 . Sub- 
sequently, Bardina et al. 12 performed a coarse-grid large-eddy simulation of homogeneous 
shear flow and Rogers, Moin and Reynolds 13 conducted fine-grid 128 x 128 x 128 direct 
simulations which further clarified its structure. 

The early experiments of Rose 5 and Champagne, et al. 6 , which were conducted for 
relatively weak shear rates and small elapsed times, seemed to indicate that the Reynolds 
stresses (and, hence, the turbulent kinetic energy) asymptoted to equilibrium values. This 
asymptotic state is consistent with the production-equals-dissipation equilibrium that 
had been hypothesized by Townsend 2 several years earlier. However, the integral length 
scales were still growing monotonically at the end of the Rose 5 and Champagne et al. 6 
experiments suggesting the strong possibility that an asymptotic state had not yet been 
reached. Subsequent physical experiments 8-10 and direct numerical simulations 11-13 , con- 
ducted for stronger shear rates and larger elapsed times, confirmed this. These physical 
and numerical experiments, along with alternative theoretical analyses (see Rogallo 11 and 
Tavoularis 14 ), have led to the following widely accepted physical picture of homogeneous 
shear flow: 
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(1) the turbulent kinetic energy k and dissipation rate e grow exponentially in time 
at the same rate, 

(2) the anisotropy tensor bij and the dimensionless turbulent time scale Sk/e reach 
equilibrium values that are relatively independent of the initial conditions and the shear 
rate S. 

It should be noted that Speziale and MacGiolla Mhuiris 15 have recently shown that 
virtually all of the commonly used two-equation turbulence models and second-order 
closures are consistent with this hypothetical picture of homogeneous shear flow. 

An alternative physical picture of homogeneous shear flow is presented in this paper 
which is consistent with the previously conducted physical and numerical experiments, 
yet at the same time is more intuitively satisfying in that it excludes the occurrence of 
unbounded energy growth. In particular, it will be shown that when vortex stretching is 
accounted for in the dissipation rate transport equation, a production-equals-dissipation 
equilibrium results in which the turbulent kinetic energy and dissipation rate eventually 
asymptote to bounded values. Illustrative calculations are presented for a k — e model 
suitably modified to account for vortex stretching. Consistent with physical and numeri- 
cal experiments, these calculations indicate an exponential time growth of the turbulent 
kinetic energy and dissipation rate for St < 30. This is the largest elapsed time considered 
in any of these previous experimental studies. However, for St > 30, vortex stretching 
eventually takes over causing the system to saturate and attain an equilibrium structure 
with bounded kinetic energy and dissipation. In the sections to follow, a detailed case is 
made to establish that this alternative equilibrium structure of homogeneous shear flow 
is a serious possibility that could have important implications for turbulence modeling 
at high Reynolds numbers. 

2. THEORETICAL BACKGROUND 


We will consider the incompressible and isothermal turbulent flow of a viscous fluid. 
The governing field equations are the Navier-Stokes and continuity equations given by 


dvi 


dvi 


dt ^ V * dxj 


dP 

dx{ 


+ 1 /V 2 Vi 



( 1 ) 

( 2 ) 
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where v,- is the velocity vector, P is the modified pressure and v is the kinematic viscosity 
of the fluid. As in the usual treatments of turbulence, the velocity and pressure will be 
decomposed into ensemble mean and fluctuating parts, respectively: 


Vi=Vi+Ui, P = P + p. 


(3) 


For homogeneous shear flow, the mean velocity-gradient tensor takes the form 


du~ 


dxj 


0 5 0 \ 
0 0 0 
0 0 0 / 


(4) 


In direct numerical simulations of homogeneous shear flow, an initially isotropic turbu- 
lence is subjected to the constant shear rate S and its time evolution is then computed. In 
laboratory experiments, an initially decaying isotropic turbulence, created downstream 
of a grid, is subjected to a uniform shear rate as it evolves spatially. The two problems 
are related, in an approximate sense, by the Galilean transformation 


X = Xo + U c t 


(5) 


where U c is a characteristic mean velocity that is typically taken to be the centerline 
mean velocity of the uniform shear [hence, dimensionless time St = 5'(s — x 0 )/U c \. 
For the remainder of this paper, we will only consider the temporally evolving version 
of homogeneous shear flow, since it is the only version of this problem that is exactly 
homogeneous. 

In any homogeneous turbulent flow, the exact transport equations for the turbulent 
kinetic energy k = \uiu~ and dissipation rate e = take the form 16 


where 


k = V-e 

(6) 

C — Ve s + 'Pey ~ 

(7) 

t> Qv\ dvl 

v = - Ti ’5v v “ = 2 " 

(8) 

„ . dui - n , du> t Doj. 

V n =2, UiUj —, *. = 2» , 

(9) 


are, respectively, the production of turbulent kinetic energy, the production of dissipation 
by mean strains, the production of dissipation by vortex stretching, and the destruction 
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of dissipation. In (8) - (9), r^- = t i;Uj is the Reynolds stress tensor and u> = V x u is the 
fluctuating vorticity vector. At this stage, we introduce the anisotropy tensor bij defined 

( 10 ) 


( t.j - § kSa) 
bi ’ = 2k 


which will be useful in the analysis to follow. 

If we non-dimensionalize the turbulent kinetic energy, dissipation rate, and time as 
follows: 

k * = T-, e* a t* = St, 
k a e 0 

where k a and e a denote initial values, then the exact transport equation (6) for k can be 
rewritten in the dimensionless form 


k* = 2 (JL _ ij b 12 k * 


( 11 ) 


since V = —T12 S in homogeneous shear flow. It should be noted that the ratio of 
production to dissipation Vjt can be related to i> 12 through the identity 


V 9 h Sk 

— = — 26 x 2 — • 

e e 


( 12 ) 


Consequently, if any two of the quantities P j e, b i2 and Sk/e are known, the remaining 
one can be computed using (12). Furthermore, since for homogeneous shear flow Vjt > 0 
and Sk/e > 0, it follows from (12) that b u < 0. Physical and numerical experiments 
have tended to indicate that b l2 and Sk/e reach equilibrium values that are relatively 
independent of the initial conditions. More specifically, these experiments suggest the 
results 

(^) « 6.0, ® -0-15, (j)^ m 1.8 (13) 

where ( . ) QO is the equilibrium value obtained in the limit as f — » oo. From (11) and (12) 
it follows that for t* 1 


where the growth rate A is given by 

k* ~ exp(At*) 

(14) 

a = 2 [(. 

».-■] 

(612)00 w 0.14. 

(15) 

Equations (13) and (14) then imply that for t* » 1 



e* ~ exp(Af*) 

(16) 
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since, if k grows exponentially and Sk/e equilibrates, then e must grow exponentially at 
the same rate. This constitutes an alternate derivation of the Tavoularis 14 asymptotic law 
of exponential growth in shear flow. The wealth of experimental data on homogeneous 
shear flow collected over the past two decades appears to be in agreement with this 
hypothetical physical picture. This is evidenced by Figure 1 containing several sets of 
measured k data versus non-dimensional time that was recently compiled by Rohr, et 
al. 10 . 


In order for this structural equilibrium — with an unbounded exponential time growth 
of k and e — to be valid, the higher order correlations in the dissipation rate transport 


equation (7) must satisfy certain consistency conditions. When non-dimensionalized, (7) 
takes the form 

/ O. . . . fit n 

e 


e* = 


2(ViU>2 


+ 




w k u> k (wjtWfc)? Sk 


\[R t - 


du>{ dufj 
2l / dxj dxj 


S w k u k 


(17) 


where Rt = k 2 jue is the turbulence Reynolds number. In deriving (17), we have made 
use of the fact that e — vuj 1 uj 1 since the turbulence is homogeneous. The first of the three 
correlations 


U)\U>2 


dui 

“Wats 


2i/ 


dwi 


WfcCUfc ’ (u>fclUfc)5 ’ W k W k 

appearing in (17) is bounded by the Schwarz inequality. Direct numerical simulations of 
Rogers 17 for homogeneous shear flow, as shown in Figure 2, indicate that the second of 
these correlations (which is proportional to the skewness Sjc in an isotropic turbulence) 
asymptotes fairly quickly to an apparent equilibrium value of 0.1. Hence, assuming the 
equilibration of ej Sk, it follows that in order to recover the exponential growth law (16) 
for f* 1, we must have the third of these correlations behave as 


dwj dui . 

v 9x i ax * ~ \[Rt ~ exp(-Af). (18) 

v 2 

While (18) may appear to be consistent with the traditionally accepted scaling laws for 
equilibrium turbulent flows 16 , it is a somewhat questionable possibility for homogeneous 
shear flow where one might expect this correlation to equilibrate. In particular, this 
term represents the ratio of the destruction of enstrophy to the enstrophy. The analo- 
gous quantity for the turbulent kinetic energy — namely, the ratio of the “destruction of 
turbulent kinetic energy” (e) to the turbulent kinetic energy (k) — does reach an equilib- 
rium value for large times. Direct numerical simulations of homogeneous shear flow do 
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not support (18), although the issue is still far from settled (see Figure 3 obtained from 
Rogers 17 , which is not suggestive of a \/Rt growth). 

We will now consider possible alternative equilibrium states for homogeneous shear 
flow. There are two such possibilities: 

(1) an alternative structural equilibrium where ( e/Sk )oo = 0, or 

(2) a production-equals-dissipation equilibrium with bounded energy states (i.e., with 
A;,*, < oo and < oo). 

The former structural equilibrium where (e/Sk)^ = 0 has been shown by Speziale 
and MacGiolla Mhuiris 15 to be primarily associated with solutions for k and e that 
undergo a power law growth with time (i.e., for St >> l,k ~ f“, e ~ t P where a > 
/3). However, these solutions are largely unstable within the context of Reynolds stress 
transport models. Furthermore, at the end of all of the previously conducted physical and 
numerical experiments on homogeneous shear flow, e/Sk either appeared to equilibrate to 
a non-zero value or continued to grow — results that are not suggestive of an equilibrium 
state where ( e/Sk ) 00 — 0. Hence, we feel that this alternative equilibrium structure is 
not a strong possibility. 

The second possibility — namely, the production-equals-dissipation equilibrium — 
appears to contradict the physical and numerical experiments which indicate that the 
turbulent kinetic energy and dissipation rate are still growing at the end of the experi- 
ments (i.e., for elapsed times St as large as 28). However, a production-equals-dissipation 
equilibrium wherein 


Voo = Coo 

(19) 

= Pcs "h P tv 

c Soo ' C V OO 

(20) 


(with bounded values for fcoo,^ and R t „) is consistent with the ensemble averaged 
Navier-Stokes equations. In the following sections, it will be shown that when the vortex 
stretching term is accounted for in the dissipation rate transport equation, then this 
production-equals-dissipation equilibrium (with an early exponential time growth of k 
and e for St < 30) becomes a serious possibility. 

3. THE DISSIPATION RATE TRANSPORT EQUATION WITH VORTEX 
STRETCHING 

Batchelor and Townsend 18 have shown that the transient behavior of the enstrophy 
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u/ 2 = u>iU) t in isotropic turbulence is governed by the equation 


dw 2 7 
dt ~ 3\/l5 


u> 3 Sk 


14 3 G 

7 =UJ 3 —- 

3i/l5 R>. 


( 21 ) 


where 
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is the skewness of the probability density function of defined with a negative sign 
to make it positive definite. In (21), R\ = u Tm3 \/u is the turbulence Reynolds number 
based on the Taylor microscale and G is a function defined by 




G = \ 4 ft v . 


(23) 


Here, u rm , is the root-mean- square of the velocity fluctuation component u, A is the 
Taylor microscale derived from the two point longitudinal velocity correlation function 
/(r), and /’ v is the fourth derivative of f(r) evaluated at r = 0. The first term on the 
right-hand side of (21) accounts for the effect of vortex stretching and is positive definite, 
while the second term is always negative and leads to the destruction of enstrophy. 

Equation (21) may be converted into a transport equation for e through use of the 
identity 

e = i/tu 2 , (24) 


which is valid for any homogeneous turbulence. Since the defining equation for the 
microscale may be conveniently written as 

, 101b 

or = 

A 2 ’ 

it follows that (21) can be converted to the equation 



_ L G t 

3\/T5 \fv 15 k ' 


(26) 


Once values are obtained for Sk and G, (26) may be solved together with the kinetic 
energy equation 

k = -e (27) 


yielding a solution for isotropic decay. 
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Batchelor and Townsend 18 showed that if G has the form 


G = y + l -RxS K , 


( 28 ) 


then the solution of (26) - (27) is compatible with their experimental data showing a 
power law decay of the kinetic energy, with an exponent of approximately one, for the 
case of moderately large values of R\. More precisely, the substitution of (28) into (26) 
yields the equation 



(29) 


which, when combined with (27), gives rise to the exact solution 


‘“•('+S" 

for isotropic decay (i.e., a power law decay where k ~ t -1 ). Apart from the specific value 
of the numerical coefficient on the right-hand side of (29), which is more commonly set 
to a value ranging from 1.83 - 1.92 to reflect more recent decay data 19 suggesting that 
k ~ f -1-1 or k ~ t -1,2 , the form of (29) is one which forms a cornerstone for the standard 
modeled e transport equation that is widely used in turbulence models. 

In the approach of Batchelor and Townsend, the choice of G given by (28) forces the 
vortex stretching term in (26) to be exactly subsumed by the action of the destruction 
of enstrophy term. This assumption, as has been noted above, guarantees compatibility 
with their isotropic decay data. However, it seems highly unlikely that two terms which 
represent distinct physical processes would exactly counterbalance each other, especially 
since the destruction of enstrophy term vanishes in the limit of zero viscosity whereas 
the vortex stretching term does not. In fact, the survival of the vortex stretching term 
in the limit of zero viscosity (i.e., the limit as R t — * oo) is crucial for the prediction of 
enstrophy blow-up — a property of solutions of the Euler equation (see Lesieur 20 ). This 
can be seen by setting v equal to zero in (21) which yields the equation 


dw 2 7 
~dT ~ 3\/T5 


3 c(°) 




K 


(30) 


where S $ is the zero-viscosity skewness. For constant > 0, (30) predicts that the 
enstrophy blows up at the critical time 


tc 


6\/T5 1 
~s® 
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where u/ 2 is the initial enstrophy. Although EDQNM supports this result, a finite-time 
enstrophy blow-up has not been seen in direct numerical simulations of the Euler equation 
(see Lesieur 20 and Pumir and Siggia 21 ). This means that either Sjp is a very small 
constant or a monotonically decreasing function of time (in the former case the enstrophy 
would blow-up at t c u> 0 whereas in the latter case it would just grow monotonically 
without bound as t — * oo). While this issue has not been fully resolved, one thing is 
clear: is not identically zero. If = 0, then (30) yields the erroneous prediction 

u) 2 — constant 


which is not supported by direct numerical simulations or theoretical analyses of the Euler 
equation (these indicate that the enstrophy grows dramatically due to vortex stretching). 
The standard modeled dissipation rate equation (29) reduces to 

dw 2 


dt 


= 0 


in the limit of zero viscosity and, hence, incorrectly predicts (in agreement with the 
4°> = 0 case) that the enstrophy is conserved. It is a simple matter to show that (28) is 
equivalent to the alternative scaling 


G ~ constant, Sk ~ R\ l 


(31) 


in the e transport equation for R\ 1. The scaling of G in (31) can be obtained by 
making the assumption that 

O . dwj dwj 

dxj dxj € 

— — ■ <x T 

(jJk^k A; 

(namely, by the assumption that the time scale of the destruction of dissipation is pro- 
portional to the time scale of the destruction of turbulent kinetic energy). While this 
scaling for G is quite reasonable, the scaling of Sk in (31) is not acceptable since it gives 
rise to a vanishing skewness in the limit of infinite turbulence Reynolds numbers (or zero 
viscosity). Hence, we will maintain the first part of equation (31) but we will not set 
to zero. A Taylor series expansion of Sk in the variable R^ 1 yields 

S K = S^ + SPR ? + ... (32) 

where S'P = \dS k / dR^ l ] R -i =0 . By substituting the first part of (31) and (32) into (26), 
we obtain the alternative modeled dissipation rate equation 

7 c(°) , -2 

t= —= j= & — C e3 — 

3VT5 sfi k 
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where (7 C2 can be taken to be a constant for R\ 1 (in deriving (33) we have made use 
of the fact that R ^ ~ R^ 1/2 and that e 3/2 /*/v = R] /2 e 2 /k). In the limit as S$ -* 0, 
the traditional modeled transport equation for e is recovered. We take C tl = 1.90, 
rather than 2,0, since this value yields a power law decay in isotropic turbulence where 
k ~ t” 1 - 1 — a value that is more in the range of the most recent experiments 19 . As 
alluded to earlier, since direct numerical simulations of the Euler equation have failed to 
yield finite time enstrophy blow-up, the zero viscosity skewness is likely to be small. 
Lagrangian based simulations of the Euler equation conducted recently by Girimaji and 
Pope 22 yielded values for S that were of the order of 10 2 . Hence, for all of the 
illustrative calculations of this paper, we will set 

= 0.01 


for simplicity (or course, S ^ could more generally be a monotonically decreasing function 
of time). 

In figure 4(a), a solution of (33) for isotropic decay is given for R to - 300 (it should be 
noted that the new e-transport equation with vortex stretching depends on R to whereas 
the more traditional model does not). It is clear that for this value of R to , vortex stretch- 
ing has little effect on the solution. In fact, for R to < 10,000 - which includes virtually 


all of the physical and numerical experiments that have been conducted on isotropic tur- 


bulence — the differences are not major. It is only for extremely high turbulence Reynolds 
numbers (i.e., for Rt a ^ 100, 000) that vortex stretching makes a major difference. For 
example, for Rt„ = 10 s , the enstrophy predicted by this new model exhibits a dramatic 
early time growth which is reminiscent of the prelude to enstrophy blow-up predicted by 
EDQNM for large values of R to (see figure 4(b) and figure VI - 4 of Lesieur 20 ). 

For anisotropic homogeneous turbulent flows, (33) can be generalized to the form 


e = n s + 


7 c(°) 

j jK a _ n — 

3v/15 vV ■>*’ 


(34) 


where V ts is the production of dissipation by mean strains given by (8). In order to 
achieve closure, a model for V Cs is needed. For simple homogeneously strained turbulent 
flows, it can be assumed that 


V es oc V 


(35) 


which, after invoking elementary dimensional analysis, yields the model 


e=C €lT V + 


e _ 

k' 


} _C — 

3s/l5 


(36) 
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where C tl is a dimensionless constant. The model for V Cs in (36) has been used in the 
k — e model of turbulence as well as in more complex second-order closures. Its success is 
largely tied to the fact that it constitutes a good approximation for plane shear flows 13,15 
— the type of flow being considered in this study. For practical calculations, C Cl can be 
taken to be 1.45 (a value obtained from equilibrium shear flows). In the next section, we 
will apply this model to homogeneous shear flow. 

4. ILLUSTRATIVE CALCULATIONS 


In order to illustrate the effect of vortex stretching on homogeneous shear flow, we 
will present the results of calculations with a k — e model for which the Reynolds stress 
tensor is modeled by 


k 2 

*ii = ~ C »~ 


dv l dvj 
9xj + dx. 


(37) 


where C p = 0.09 is a dimensionless constant. While the k — e model is somewhat 
simplistic since it is based on an eddy viscosity, it was recently shown by Speziale and 
MacGiolla Mhuiris 15 that this model is topologically equivalent to the more complex 
second-order closure models for homogeneous shear flow (the deficiencies in the k — e 
model only become pronounced when there are combinations of shear and rotation or 
multi-dimensional strains). Hence, (37) will suffice to illustrate the qualitative changes 
induced when the effect of vortex stretching on the dissipation rate is accounted for. 
Equation (37) will be solved in conjunction with equations (6) and (36). The standard 
k — e model is recovered in the limit as 4 0) - 0. 

For homogeneous shear flow, the k — e model with vortex stretching yields the trans- 
port equations 


k 2 

k = C u —S 2 - e 


(38) 


and 


e = 


C Cl C tl kS 2 + 


7 * e 2 

3v/l5 \/v “ £J k 


(39) 


which are obtained by substituting (4) into (6), (36) and (37). 

For all non-zero values of S$ it is a simple matter to show that the solution to (38) 
and (39) converges to an equilibrium state, with bounded energies. The equilibrium 
values may be found by setting the right-hand-sides of (38)-(39) to zero, yielding the 
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results 


and 


k 0 o 135 — C ei ) 2 Sk 0 


49 


c(°r v> 

Jjc U tc 


^OO 

Co 


135C„(C7 Q -C tl ) 2 


49 


7 (o)' 

>K 


Ri 


' Sko ' 
> € 0 


(40) 


(41) 


These relations clearly indicate that A^/A^ and e^/e,, have a SpP dependence so that 
the standard k — e model prediction of an unbounded growth of k and e is easily recovered 
in the limit as S$ — > 0. 

Using (40) and (41) the following additional equilibrium values are also obtained for 
this k — e model with vortex stretching: 


and 


- - In' 

\ 6 J oo 

(42) 

fir). - v® 

(43) 

l35(C e2 -C £l ) 2 

49 e(°) 2 

(44) 




These results (which are independent of the initial conditions) differ from the values of 
(Sife/e)*, = y/a/C^, (-m/k) oo = yfaC^ t and R too = oo obtained from the standard 
k - e model where a = (C ej - l)/((7 ei - 1) « 2. The mechanism by which the presence 
of the vortex stretching term has the effect of creating bounded long time solutions lies 
in its enhancement of the growth rate of e. Evidently, this increase in the growth rate 
of e is accompanied by a simultaneous reduction in the production of k , thus forcing 
a production-equals-dissipation equilibrium. In alternative terms, vortex stretching — 
which becomes more pronounced at high turbulence Reynolds numbers since it scales as 
— eventually causes homogeneous shear flow to undergo a saturation to an equilib- 
rium state with bounded component energies (the values of which are set by the shear 
rate and the viscosity). 

We will now show that the k — e model with vortex stretching yields temporal evo- 
lutions of the turbulence fields for St < 30 that are in good qualitative agreement with 
previously conducted physical and numerical experiments. When non-dimensionalized, 
(38) - (39) take the form 

Sk Q k t a , 


k* = Co 


Sk 0 


(45) 
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(46) 


f = C„C„— *• + 




e„ 3VI6 K Sk„v-- — SKk • 

where again we have (7 M = 0.09, (7 Cl = 1.45, (7 €2 = 1.90 and 5^ = 0.01 (the standard 
A? — e model is obtained by setting SfP = 0). The initial conditions, which correspond 
to an isotropic turbulence, are taken to be e 0 /Sk 0 = 0.296 and R to = 300. These are 
the approximate initial conditions of the large-eddy simulations of Bardina et al. 12 which 
will allow us to make some direct comparisons between the model and the simulations. 

Figures 5(a) and 6(a) display the short time solutions for k* and e* compared to the 
large eddy simulation data. The solutions with non-zero vortex stretching are seen to 
display short term exponential growth in a manner very similar to that in the standard 
k — e closure. The effect of the vortex stretching term is to reduce the growth rate of 
k* and e*, though initially there is a slight increase in the magnitude of e*. A view of 
these solutions over a much longer time interval, as displayed in Figures 5(b) and 6(b), 
reveals the dramatic effect that the vortex stretching term ultimately has on the long 
term growth of k* and e*. It is seen that with the vortex stretching effect included, the 
initial exponential growth rates are eventually suppressed, so that by St ^ 40, k* and e* 
asymptote to bounded equilibrium values. 

Figures 7 and 8 show the short and long time behavior of the dimensionless ratios Sk/e 
and —TFo/k for the solutions obtained with and without vortex stretching. The curves in 
Figures 7(a) and 8(a) give the impression that an equilibrium state for these quantities 
may have been achieved by the time St 10. However, the long time solutions reveal 
that, in the case where vortex stretching effects are included, only a local maximum is 
reached — further developments must occur before a true equilibrium is achieved. This 
suggests that the apparent convergence of quantities such as Sk/e seen in numerical and 
experimental studies, may not signify that a final equilibrium state has resulted (i.e., it 
may only be a local maximum). The behavior of Sk/e shown in Figure 7(b) with vortex 
stretching present is remarkably consistent with the direct simulations of Rogers et al. 13 
shown in Figure 9 (i.e., the tails in the computations of Rogers et al. 13 may indicate that 
an equilibrium state has not been reached). 

Figure 10 provides a plot of the long time behavior of the computed turbulence 
Reynolds number. It achieves an equilibrium value of approximately 5600, which is more 
than eighteen times its initial value. The most significant effect of the vortex stretching 
term on lb* and e*, as seen in Figures 5(b) and 6(b), occurs for Rt > 4000. This confirms 
the belief that vortex stretching is mostly a phenomenon associated with high turbulence 


13 


Reynolds numbers. Finally, some indication of the sensitivity of the computed solutions 
to the skewness 5^ appearing in the vortex stretching term is shown in Figure 11. This 
contains the time evolution of k* for a range of values of as expected, k ^ increases 
with decreasing values of Sg\ Figure 12 shows the effect on k* of a change in the 
initial values of e Q /Sk. As would be expected on physical grounds, an increase in the 
dimensionless shear rate leads to a higher equilibrium value for the turbulent kinetic 
energy. 

5. CONCLUSIONS 

An alternative view concerning the equilibrium structure of homogeneous turbulent 
shear flow has been presented based on the inclusion of the effect of vortex stretching. It 
was shown that the presence of just a small amount of vortex stretching can ultimately 
lead to a production-equals-dissipation equilibrium with bounded energy states. For 
elapsed times St < 30 — which includes the largest values of St considered in any of the 
previously conducted physical and numerical experiments — the introduction of a small 
vortex stretching term into the standard modeled dissipation rate transport equation still 
yields an exponentially growing kinetic energy and dissipation for early times. However, 
since this vortex stretching term scales as \/Rty it eventually becomes dominant, causing 
a saturation of the system to an equilibrium state with bounded turbulent kinetic energy 
and dissipation. Although this alternative physical picture of homogeneous shear flow is 
contrary to the commonly accepted asymptotic laws, for which an unbounded exponential 
time growth of k and e is postulated, it is a real possibility that needs to be seriously 
considered in the future. New physical and numerical experiments, for larger elapsed 
times St } could shed more light on the issue. The recent experiments of Tavoularis 
and Karnik 9 , which were conducted up to St = 28, did show some tendency of the 
integral length scales to level off — a feature which, if more solidly established, would 
be supportive of the existence of a production-equals-dissipation equilibrium. However, 
it is possible for the kinetic energy and dissipation rate to equilibrate while the integral 
length scales grow without bound (perhaps this is the more likely possibility since the 
flow field is infinite). 

Finally, some comments should be made concerning the implications that the results 
of this paper have on turbulence modeling. Since all of the commonly used two-equation 
models and second-order closures based on the turbulent dissipation rate neglect vortex 
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stretching, they predict an unbounded exponential time growth of k and e in homo- 
geneous shear flow. This type of behavior has been shown to cause problems in the 
calculation of certain inhomogeneous turbulent flows. The singularity in plane stagna- 
tion point turbulent flow represents a prime example 23 . Hence, if it is ultimately justified, 
the alteration of turbulence models to yield a production-equals-dissipation equilibrium 
in homogeneous shear flow via vortex stretching could make their behavior more robust 
in other turbulent flows. A full resolution of this issue is only likely to come from a rig- 
orous mathematical proof based on some appropriate energy method. For the meantime, 
however, the results of this study clearly establish the need to re-examine this issue. 
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Figure 4. Comparison of solutions for isotropic decay: 5^ = 0.01; - - - - standard 

k — e model, (a) Time evolution of the turbulent kinetic energy for R te = 300, and (b) 
time evolution of the enstrophy for Rt 0 = 10 6 . 
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Figure 4. Comparison of solutions for isotropic decay: S # = 0.01; - - - - standard 

k — e model, (a) Time evolution of the turbulent kinetic energy for R te = 300, and (b) 
time evolution of the enstrophy for R to = 10 6 . 
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Figure 5. Comparison of the model predictions with the large-eddy simulations of Bar- 
dina, et al. 12 : o k * from large eddy simulation 12 ; k* obtained with SjP = 0.01; - - - 

- k* obtained from the standard k — e model, (a) Short time solution, and (b) long time 
solution. 
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Figure 5. Comparison of the model predictions with the large-eddy simulations of Bar- 
dina, et al. 12 : o k* from large eddy simulation 12 ; k* obtained with 5^ = 0.01; - - - 

- k* obtained from the standard k — e model, (a) Short time solution, and (b) long time 
solution. 
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(a) 



Figure 6. Comparison of the model predictions with the large-eddy simulations of Bar- 
dina, et al. 12 : o e* from large eddy simulation 12 ; e* obtained with SjP = 0.01; - - - 

- e* obtained from the standard k — e model, (a) Short time solution, and (b) long time 
solution. 
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Figure 6. Comparison of the model predictions with the large-eddy simulations of Bar- 
dina, et al. 12 : o e* from large eddy simulation 12 ; e* obtained with SjP = 0.01; 

- c* obtained from the standard k — e model, (a) Short time solution, and (b) long time 
solution. 






Figure 7. Comparison of the model predictions for Sk/e: 5 

k — e model, (a) Short time solution, and (b) long time solutio 
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Figure 10. Time evolution of the turbulence Reynolds number Rt predicted by the k — 
model with vortex stretching. 
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Figure 11. Sensitivity of the model predictions for the turbulent kinetic energy to S$: (a) 
S { JP = 0.0133; (b) Sp = 0.01; (c) SP = 0.0066; (d) S { P = 0.005 and (e) sP = 0.0033. 


32 








NASA 

Natooai Ae^onauics and 
Space Aorntnisiraion 


Report Documentation Page 


1. Report No. 

NASA CR- 187439 
ICASE Report No. 90-66 


2. Government Accession No. 


3. Recipient's Catalog No. 


4. Title and Subtitle 

BOUNDED ENERGY STATES IN HOMOGENEOUS TURBULENT 
SHEAR FLOW — AN ALTERNATIVE VIEW 


5. Report Date 

October 1990 


6. Performing Organization Code 


7. Author(s) 

Peter S. Bernard 
Charles G. Speziale 


8. Performing Organization Report No. 

90-66 


9. Performing Organization Name and Address 

Institute for Computer Applications in Science 
and Engineering 

Mail Stop 132C, NASA Langley Research Center 
Hampton, VA 23665-5225 


10. Work Unit No. 


505-90-21-01 


11. Contract or Grant No. 

NAS1-18605 


12. Sponsoring Agency Name and Address 
National Aeronautics and Space Administration 
Langley Research Center 
Hampton, VA 23665-5225 


13. Type of Report and Period Covered 

Contractor Report 


14. Sponsoring Agency Code 


15. Supplementary Notes 
Langley Technical Monitor: 
Richard W. Barnwell 


To be submitted to Physics of Fluids 
A. 


Final Report 

16, Abstract 


The equilibrium structure of homogeneous turbulent shear flow is investigated from a theoreti- 
cal. standpoint. Existing turbulence models, in apparent agreement with physical and numerical 
experiments, predict an unbounded exponential time growth of the turbulent kinetic energy and dis- 
sipation rate; only the anisotropy tensor and turbulent time scale reach a structural equilibrium. 
Tt is shown that if vortex stretching is accounted for in the dissipation rate transport equation, 
then there can exist equilibrium solutions, with bounded energy states, where the turbulence pro- 
duction is balanced by its dissipation. Illustrative calculations are presented for a k-e model 
modified to account for vortex stretching. The calculations indicate an initial exponential time 
growth of the turbulent kinetic energy and dissipation rate for elapsed times that are as large as 
those considered in any of the previously conducted physical or numerical experiments on homoge- 
neous shear flow. However, vortex stretching eventually takes over and forces a production-equals- 
dissipation equilibrium with bounded energy states. The validity of this result is further sup- 
ported by an independent theoretical argument. It Is concluded that the generally accepted struc- 
tural equilbrium for homogeneous shear flow with unbounded component energies is in need of re-ex- 
amination . 


17. Key Words (Suggested by Author(s)} 
turbulence; homogeneous shear flow; 
vortex stretching; K-e model 


18. Distribution Statement 

34 - Fluid Mechanics and Heat Transfer 



Unclassified - 

- Unlimited 


19. Security Classif. (of this report) 
Unclassified 

20. Security Classif. (of this page) 
Unclassified 

.. ... 

21. No. of pages 

35 

22, Price 

A0 3 


NASA FORM 1626 OCT 86 


NASA-Langley, 1990 



